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A commonly adopted relational account of time evolution in generally-covariant systems, and more 
specifically in quantum cosmology, is argued to be unsatisfactory, insofar as it describes evolution 
relative to observed readings of a clock that does not exist as a bona fide observable object. A 
modified strategy is proposed, in which evolution relative to the proper time that elapses along the 
worldline of a specific observer can be described through the introduction of a 'test clock', regarded as 
internal to, and hence unobservable by, that observer. This strategy is worked out in detail in the case 
of a homogeneous cosmology, in the context of both a conventional Schrodinger quantization scheme, 
and a 'polymer' quantization scheme of the kind inspired by loop quantum gravity. Particular 
attention is given to limitations placed on the observability of time evolution by the requirement 
that a test clock should contribute only a negligible energy to the Hamiltonian constraint. It is 
found that suitable compromises are available, in which the clock energy is reasonably small, while 
Dirac observables are reasonably sharply defined. 

PACS numbers: 98.80.Qc, 04.20.Cv, 04.60.Ds 



I. INTRODUCTION 



o . 

Taken at face value, the canonical formulation of general relativity seems to entail that time evolution is equivalent 
$H ' to a gauge transformation, and should therefore be physically unobservable. Not surprisingly, this 'problem of time', 
along with other related interpretational issues, has attracted considerable attention over many years. The situation in 
the early 1990s was comprehensively documented in 0,0; a recent survey is presented in Q, and textbook discussions 
{SJ . can be found in Everyday observations suggest that things actually do change with time, and most investigators 

£> ' have concluded, in one way or another, that this fact can be accomodated only by reference to physical clocks. Many 
different schemes for implementing this general idea have been proposed, but the notion of 'relational', 'emergent' 
or 'internal' time, originating in the work of Rovelli [H-IH has been quite widely adopted, especially in the context 
of quantum cosmology, where many explicit calculations have become possible in recent years (see, e.g. (9rlll| for 
reviews). Roughly speaking, the relational-time approach involves identifying, within the model considered, some 
quantity that is to serve as a clock (in simple cosmological models, this is typically a scalar field) and describing the 
evolution of other quantities relative to putative values of this clock variable, which thus serves as an 'internal time'. 
In this paper, we first argue that such accounts of time evolution, while mathematically sound, are hard to interpret 
in a consistent manner, and cannot be the whole story. We then propose a modified account, which we of course 
believe to be less unsatisfactory. 

For concreteness, we discuss in section [TTJ a particular implementation of the notion of internal time given in [121414} 
(an amended version of the last of these papers appears in [l5j . and a rather different implementation, based on 
effective semiclassical dynamics is proposed in fl6j ) and explain in some detail why we believe it to be deficient. 1 In 
brief, this is because (i) the physical clock does not appear as an observable object in the final theory, so its putative 
readings cannot be interpreted as values obtained by observation and have, indeed, no clear meaning; and (ii) no 
account is available of the passage of time as experienced by real observers. 

In sections Mil and l lVI we describe, in the context of two different quantization schemes, how these difficulties might 
be overcome, through the introduction of a 'test clock' associated with some specific observer (an idea proposed some 
time ago in |17|). This makes it possible to describe time evolution from the point of view of a particular observer, with 
respect to the proper time that elapses along that observer's worldline. A test clock should serve to reveal the time 
evolution of the model universe while contributing negligibly to its energy content, and this restricts the resolution 
with which the values of time-dependent observables can be determined. For both the conventional 'Schrodinger' 
quantization scheme of section Mil and the loop-quantum-gravity-inspired scheme of section IIVI we investigate this 



'Electronic address: i.d.lawrie@leeds.ac.uk 

1 We emphasize that this discussion is not intended as a criticism of the work presented in Il2l4l5ll . On the contrary, these papers describe 
an especially clear and thorough implementation of a view of time evolution that is now sufficiently widespread to have something of 
an 'official' status. We do, of course, wish to suggest that this 'official' view is not wholly satisfactory. 
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restriction, and find that suitable compromises are possible. This proposal can be seen as a variant of the general 
idea of relational time, and we discuss its relation to other versions in section [V] 



II. RELATIONAL TIME IN QUANTUM COSMOLOGY 

The homogeneous, spatially flat cosmological model with no cosmological constant studied in fl2T - fl5l ] can be specified 
classically by the action 



S = / dsN(s) 



1 (d s v) 2 1 2 
v{d s (j)Y 



24ttG v 2 



(2.1) 



Here, v(s) = a 3 (s), where a(s) is the usual Robertson- Walker scale factor, N(s) is a positive, but otherwise arbitrary 
lapse function and s a correspondingly arbitrary time coordinate, while the matter content is represented by the 
massless scalar field <fr(s). We take this action to refer to a fiducial cell of coordinate volume J d 3 x — I, and v, with 
dimensions (length) 3 to be the physical volume of this cell. In SI units, G = Gnc~ 2 , where Gn is the usual Newton 
constant. The momenta conjugate to v and <fi are 

I v 



( 2 - 2 ) 
1 27rG v 

P4> = W>> (2-3) 
where the overdot denotes differentiation with respect to the proper time 

t(a)= [ S N(s')ds', (2.4) 



so v — N 1 dv/ds, etc. Variation of the action (|2.II) with respect to N yields the constraint Ho = 0, where 

H = -6nGvp 2 v + \v~ X v\, (2.5) 

is the generator of translations in t. 

The problem of describing time evolution arises, as is well known, from the fact that bona fide gauge-invariant 
observables (Dirac observables) must commute with the constraint, and must therefore be constants of the motion. 
The solution to this problem is often taken to involve identifying a variable within the theory which can serve as a 
physical clock, and thus provide a notion of evolution with respect to an 'internal' or 'emergent' time. In [T^ - fT5j j . this 
relational picture is implemented in the quantum theory by casting the constraint in the 'deparametrized' form 



'* m ^: v) = -\l®(v, Pv M<P, «), (2.6) 

where V0 is the square root of a suitably defined operator corresponding to the classical expression 12itGv 2 p1. The 
precise form of this operator depends on the particular quantization scheme adopted. This is formally similar to a 
non-relativistic Schrodinger equation 

at 

and appears to provide a notion of evolution with respect to an internal time represented by the scalar field <j>. (A 
much earlier study of quantum cosmology, also using a massless scalar field as an internal time, was presented in [18j.) 
In particular, an operator such as 



V(fo) := cxp (-iV^O - M/h) wexp (iVI^ - O )A 



(where v acts by multiplication on ^((j), v)) is, for any fixed value of the parameter </>o, a gauge-invariant Dirac 
observable: if ty(<j),v) is a solution to the constraint equation (|2.6j) . then V{(po)^{4'^ v ) is another solution. In the 
language introduced by Rovelli, V^o) is an 'evolving constant of the motion' @, providing a 1-parameter family 
of 'complete observables' Q labeled by <f>$. For the model considered here, the classical solution for <fi is always a 
monotonic function of t, so it is tempting to interpret the constraint (|2.6p as effectively describing evolution with time, 
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in such a way that V(<t>o) represents "the volume at the time when the scalar held has the value <f>o" . (Restrictions 
on the choice of variables that might serve as 'internal time' are discussed in [Hj].) For the reasons we are about to 
present, we think that this interpretation has significant limitations, and in subsequent sections we will suggest how 
some of them might be overcome. 

The limitations we have in mind are indicated by the following interrelated observations: 

1. Despite their formal similarity, the constraint equation (12.6P and the Schrodinger equation (|2.7[) do not mean 
the same thing. In the case of a non-relativistic particle, the wavefunction ^(t, x) is, for each fixed value of t, 
an element of the physical Hilbert space "H p hys = L 2 (K, dx) corresponding to a possible instantaneous state of 
the particle. A solution of (12.71) yields a sequence of such states, labeled by the external time parameter t, and 
in that apparently straightforward sense describes the time evolution of the state of the 1-particle system. By 
contrast, the physical Hilbert space of the cosmological model is a space of solutions of the constraint equation 
(|2.6j) : a solution of this equation specifies not a sequence of possible states, but a single state characterized by 
a certain correlation between 'partial observables' v and <j>. 

2. Rovelli Q defines a 'partial observable' as a quantity for which a measurement procedure can be specified, in 
contrast to a 'complete observable', whose value can be predicted by theory. He appears to take the view that 
the time parameter t in (|2.7[) is, in this sense, a partial observable, but we disagree. The parameter t does not 
refer to the reading of any physical clock. It is an external parameter, more like Newton's 'absolute, true and 
mathematical time' or, as described by Unruh and Wald [I(J[2lj], a 'heraclitian time', which 'sets the conditions' 
for a measurement to be made. To be sure, the times recorded in a laboratory notebook during the course of 
an experiment intended to test the validity of (|2.7[) will refer to the readings of some physical clock. But then, 
according to standard quantum mechanics, the state of the combined system of a particle (position x) and clock 
(pointer reading T, say) is described by a wavefunction ^(t,x,T), governed by its own Schrodinger equation. 
Under suitable conditions, a sequence of observed values of T may closely approximate the corresponding values 
of t at which the observations were made, and one might derive an approximate version of (|2.7[) in which t is 
replaced with T . However, this description is necessarily approximate, and there are well known restrictions on 
the ability of a quantum- mechanical clock to furnish a reliable measure of t (see, e.g. [2l|-[24j]). 2 Consequently, 
if the scalar field <f> in (|2.6I) is taken as analogous to the reading T of a laboratory clock, we should expect that 
equation to be only approximately valid. That is not so, however: (|2.6p is an exact constraint equation, not an 
approximate evolution equation. 

3. If the constraint equation (|2.6[) is to be regarded as expressing evolution with respect to an internal 'time' </>, 
how is its solution, ^>(4>, v) to be interpreted? A statement to the effect that \^(<f>, w)| 2 dw is the probability of 
finding the volume to have a value between v and v + dv at 'time' <f> (or a similar statement that replaces dv with 
a more appropriate measure if necessary) as in standard quantum mechanics will not do, because the scalar field 
is a physical quantity that has no definite value until it is measured; <f> is not an external parameter with the 
heraclitian property of 'setting the conditions' for a measurement of v. For the same reason, the 'Heisenberg- 
picture' operator (|2.8p cannot be construed as representing the volume at 'time' <pQ. Nor will it be possible to 
interpret \^(4>, v)\ 2 in terms of a joint probability for obtaining the pair of values </> and v from simultaneous 
measurements of the scalar field and the volume. The reason is that there are not enough physically meaningful 
quantities available to be measured. Classically, the 4-dimensional kinematical phase space, with coordinates 
(v, p v , <j>, p$) is reduced by the constraint to a 2-dimensional physical phase space of distinct gauge orbits, with 
a single pair of conjugate coordinates. Quantum-mechanically, this means that only one independent quantity 
is available to be measured. Equivalently, one cannot define two independent operators v and <fi acting in the 
physical Hilbert space; there is only one physical degree of freedom, corresponding, for example, to the 'complete 
observable' V(0). 3 



2 More precisely, a wavefunction that realizes an exact correlation between T and t must have the form 9(t, x,T) = 8(T — t)tp(T,x); at 
each instant t, it is an eigenfunction of the pointer operator T that acts by multiplication. The Schrodinger equation admits solutions 
of this form only if the clock Hamiltonian is its conjugate momentum, Hj> = pj> [23| which, being unbounded below, is physically 
unrealistic. Moreover, this wavefunction is not in the physical Hilbert space L? (M 2 , drrdT) , so even in this idealized case, the possibility 
of identifying values of t with the results of measurements made on the clock is doubtful. No real clock will exist perpetually in a 
sequence of eigenstates of T, so a Schrodinger equation of the form icVi/>(T, x) = H(x, d x )ip{T, x) gives at best an approximate, effective 
description of the correlations exhibited by sequences of measured values x and T, ignoring, for example, the loss of unitarity resulting 
from the repeated measurements needed to obtain these values |24| . 

3 The same point can be phrased in terms of probability measures: if '!'(</>, v)\ 2 is to be interpreted as a joint probability density, it must 
be normalized with respect to a probability measure d/i(<p, v) on a 2-dimensional configuration space. However, the physical Hilbert 
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It can be argued [25j that even in non-relativistic quantum mechanics the external time t is irrelevant to physics. 
After all, a student who investigates the motion of a pendulum has only a set of recorded position measurements 
and stopwatch readings to work with in any subsequent analysis. For generally-covariant systems, the possibility 
of describing physics entirely in terms of correlations between Dirac observables has been studied in some detail in 
[26| . From a purely operational point of view, it is no doubt true that substantive physics deals only with correlations 
between measured quantities, but this seems to offer an impoverished account of the world as it is actually experienced. 
A student who forgets to bring a stopwatch to the lab is not thereby prevented from seeing a pendulum swing; and 
a statement such as "the stopwatch read 5s some three seconds after it read 2s" seems not to be entirely vacuous, 
regardless of whether it is checked for accuracy with the use of a further clock. Moreover, while the notion of time 
may be scarcely less nebulous in Newtonian mechanics than it notoriously appeared to Augustine [13], classical general 
relativity provides a concrete meaning for t as the proper time that elapses along an observer's worldline, even though 
it affords no experimental procedure for determining the actual values of t. That is, it seems meaningful to regard the 
readings of a physical clock as supplying an estimate (more or less reliable according to the quality of manufacture) of 
the intervals of geometrical proper time that elapse along its worldline, even though we have no experimental means 
of checking that this is actually so. 

To summarize, the interpretation of the wavefunction ^{<fi, v) that satisfies the constraint equation (|2.6p presents 
two related difficulties. On the one hand, there is no 'heraclitian' time variable that would allow us to make sense of 
the theory in terms of our familiar sense of the passage of time. The volume v and scalar field (f> are physical quantities 
which, in principle, can be measured. (That is, they are apparently the same sort of thing as electromagnetic fields, 
and only the practical difficulty of constructing apparatus that couples to them stands in the way of making such 
measurements.) It ought to be possible to formulate questions such as "given that I have just determined the volume 
and scalar field to be vq and cj>o, what is the probability that I will find them to be v\ and <fii in an hour's time?", 
but the theory as it stands does not admit such questions. On the other hand, if we rule that questions of this kind 
are inadmissible, and confine ourselves to studying correlations between measured values of v and <fi, we find that this 
cannot be done either, because there are not enough Dirac observables to be measured. 

Of course, a more comprehensive theory will afford more Dirac observables to be correlated, but that is not of much 
help. For example, a theory with an extra scalar field deals, apparently, with three measurable quantities, v, <p and 
tp, say, but yields only two independent Dirac observables, say V(tpo) and <!>(-0o), constructed along the lines of (|2.8|) . 
We still face the problem that the number i/jq cannot be interpreted as the result of a measurement of ip, because 
there is no corresponding Dirac observable available to be measured. Clearly, the same will apply to theories with 
more than one geometrical variable, such as the Bianchi I model studied i n 12811 . 



For the simple model studied here, the wavefunction ^{(j), v) computed in [12l.ll3| looks as if it describes a correlation 
between v and 0, and seems to be peaked along a classical trajectory (<fi{t),v(t)). The problem is that this appearance 
is at variance with the number of observables to hand (or, equivalently, with the dimension of the configuration space 
on which the wavefunction is defined — see footnote [31) 

In the remainder of this paper, we describe a possible solution to these difficulties, based on a view of time evolution 
proposed in [13] ■ It is useful to suppose that, as in standard quantum mechanics, the wavefunction ^(<j), v) describes 
the state of a system from the point of view of an observer external to the system itself. In our case, the hypothetical 
observer is external to the entire model universe, and any sense that this observer might have of a 'passage of time' 
is quite separate from what passes inside the universe. The external observer has the possibility of determining the 
value of only a single Dirac observable; this one value completely specifies the state of the universe-which means its 
entire history. The difficulties identified above do not concern the information accessible to the external observer, 
which is delivered by the wavefunction ^{4>, v) according to the usual quantum- mechanical rules, but rather what 
might be observed by an observer internal to the universe. To assess the latter, it is necessary to include in our model 
a description of relevant features of the physical system that does the internal observing: at a minimum, the clock 
from which this system gains its sense of time. By doing this, as we shall illustrate, it is possible to alleviate both 
of the difficulties. On the one hand, we introduce a genuine 'heraclitian' time r, corresponding to the arc-length 
of the observer's worldline, not to readings of a physical clock. Time evolution with respect to r is described by a 
standard Schrodinger equation, precisely analogous to (|2.7p . though in practice we will deal with the corresponding 
Heisenberg-picture operators V(t) and ^(r). (Unruh and Wald |2l| argue that the problems of interpretation are not 
alleviated by the introduction of observers, but the role they envisage for these observers is different from the one used 
here.) On the other hand, the lack of a Dirac observable corresponding to the reading of the physical clock can be 
understood by considering that this clock is internal to the observing apparatus, and thus in principle unobservable 




space to which a solution of 1 12.611 belongs is something like L 2 (R, d/i(v)), the details depending on the quantization scheme, with a 
^-independent inner product. 
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by the observer from whose point of view the evolution is described. 

Clearly, the idea of incorporating an observer's clock into our model is in some way akin to the notion of a material 
reference frame, which has been widely studied (see, for example [H 13. PR [2914351 ] ) . We think that the implementation 
of this general idea described below differs in important respects from others to be found in the literature, and will 
return to this point in section IVl 



III. CONVENTIONAL QUANTIZATION 



As a rule, cosmologists do not find it necessary to include the energy content of their observing apparatus explicitly in 
the Friedmann equation. Correspondingly, we seek to modify the cosmological model defined by (|2.1j) by introducing 
a 'test clock', which will serve to reveal the time evolution of the volume and scalar field, while disturbing the 
Hamiltonian constraint to a negligible extent. In this section, we consider the quantum theory of such a model, 
using a conventional quantization scheme similar, though not identical, to the 'Wheeler-de-Witt' scheme described 
in [12l4l5| , within which most of the analysis can be achieved exactly and explicitly. In section IIV1 we will consider 
a 'polymer' quantization scheme of the type employed in loop quantum cosmology, where we will need to resort to 
approximations. In either case, we may expect (naively, on the basis of an 'energy-time uncertainty relation') that 
restricting the energy of the clock should place some limit on the resolution with which the value of a time-dependent 
Dirac observable can be determined, and we will investigate this issue in some detail. 

A. Homogeneous cosmology with a test clock 

Classically, consider a small clock, whose internal workings are described by a Lagrangian £(X), localized on the 
worldline 4 x fJ '(X), parametrized by its arc length A. Its contribution to the action has, in general, the form 



clock 



d^vW J dX£(X)6(x,x(X)) (3.1) 
dX£(X) (3.2) 



where S is a covariant delta function, with the property J d 4 xy / '\g(x)\5(x, x')f(x) = f(x'). In principle, the coordinates 
x fJ '(X) are dynamical variables, but in order to construct a simple cosmological model, we remove these degrees of 
freedom, along with most of those in g^, by taking the world line to be that of a comoving observer in a FRW 
universe. Then the arc length A coincides with t in (|2.4I) and the total action becomes 

(3.3) 



S = J dsN(s) 
The Hamiltonian constraint now reads 



1 (dsv) 2 , 1„, Q . a 



24ttG v 2 



v(d s <f>) 2 + N( S ) 2 £(t(s)) 



H + h= -6nGvp 2 v + \v- l pl + h = Q (3.4) 

where h is the Legendre transform of £, and is intended to be very small compared with the matter term in H$, in the 
same sense that a silicon chip on board WMAP is small compared with the total energy content of the visible universe. 
Since the action (|3.3p is supposed to apply to a homogeneous universe, it might be more consistent to regard the 
last term as arising from a space-filling congruence of clocks. However, if h is small enough, then it seems very likely 
that, say, a spherically symmetric cosmology, with the clock definitely localized on a single world line, but departing 
only to this tiny extent from a genuinely homogeneous one, would be described by the action Q3.3P with negligible 
error. What is essential for our purpose is that, because the clock (or each clock) is localized on a single worldline, h 
is independent of the cosmological variables (v,p v , (j>,p$). Having obtained the constraint, we will dispense with the 
arbitrary coordinate s, and deal only with the proper time t or, equivalently, set N(s) — 1 so that s and t coincide. 



4 A more rigorous derivation than we attempt here would, amongst other refinements, consider the clock to be localized in a region 
somewhat larger than its Schwarzschild radius, but we assume that such refinements are inessential to the issue of time evolution that 
is our main concern. 
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Up to the choice of an origin, t has a clear physical interpretation as the geometrical proper time that elapses along 
the comoving observer's worldline. 

We will not be specific about the microscopic constitution of the clock. Suppose that h depends on several micro- 
scopic phase-space variables, which we denote collectively by x, and that x{x,t) is the phase-space trajectory (the 
solution of dtx — {x, h}) that passes through x at, say, t = 0. The reading of the clock is some function r(x), and 
we define the function to(x) as the solution of 

r{x{x,t )) = 0. (3.5) 

(The clock is fit for purpose only if this equation has a unique solution.) 

For a given state of the clock, the quantity r = t — to has a classical interpretation as the interval of proper time 
that has elapsed since the clock read 0. This fiducial event provides a physical meaning for the origin r = 0. By 
contrast, the arc-length parameter t is defined by (|2.4j) only up to an arbitrary constant, because a change 5N(s) in 
the undetermined lapse function changes t by an additive constant. Fixing an origin for r is the only essential role of 
the physical clock. It is not hard to show that 

{t ,h} = -l (3.6) 

and this is the only property of the clock that will matter. 

Finally, let w denote collectively the phase-space variables that appear in Ho (they are (v,p v , 4>,p^) in the present 
example, but the following result is general). For some quantity f(w), let f(w,t) be the solution of 

d t f = {f,H } (3.7) 

with the initial condition f(w,0) = f(w). Then 

F(r):=f(w,t +T) (3.8) 

is classically the value of / at a proper time r after the clock read 0. Since h and to are independent of w, and the 
vector field {• , h} is a linear differential operator, we easily find 

{F(r),h} = d T F( T ){t , h} = ~{F(t), H q }, (3.9) 

so for each / and r, F(t) commutes with the constraint H = Ho + h. Thus, F(t) gives a 1-parameter family of Dirac 
observables 5 , and obeys the equation of motion 

d T F(r) = {F(t),H }. (3.10) 

For the model at hand, it proves convenient to define y — vp v . Then the Hamilton equations are 

d t v = -12irGy (3.11) 

d t y = H (3.12) 

d t 4> = u -1 P0 (3.13) 

dtp* = (3.14) 

and they have the solutions 

v(t) = v - \2irGyt - 6nGH t 2 (3.15) 

y{t) = y + Hot (3.16) 

4>{t) = 0+_^ln(^^) (3.17) 
V127rG \v~a + tj 



5 The construction of these observables is somewhat similar to the construction of 'evolving constants of the motion' proposed by Rovelli 
0, as mentioned above, and elaborated by others, but it is not the same construction. We think the difference is important, and discuss 
it in detail in section FVl 



where a± = 67rGy ± \Z3nGprf, arise from the factorization v(t) =v 1 (v — a + t)(v — a_t). We can now define a collection 
of basic Dirac observables 

V = v - l2irGyt - QnGH tl (3.18) 

Y = y + H t (3.19) 

* = J> + ^J=ln( V -^-) (3.20) 
sfVhvG \v-a+t J 

- P0. (3.21) 
Of course, the form of the Hamiltonian is preserved 

= -QttGV^Y 2 + ^V~ l Pl, (3.22) 

the gauge-invariant variables V and Y inherit the Poisson-bracket relations satisfied by v and y, namely 

{V,Y} = V, {V,H } = -12ttGY, {Y,H } = H , (3.23) 

and the equations of motion (|3.10p have the obvious solutions 

V(t) = V - YI-kGYt - 6ttGH t 2 (3.24) 

Y(t) = Y + Hqt (3.25) 

$(r) = $ + _L= m f^4^ ( 3 - 26 ) 



Vl2nG \V -A 

P$(r) = Pi, (3.27) 



with A± = QitGY ± V3^GP* 



B. Quantum Dirac observables 



Formally, it is easy to construct quantum-mechanical versions of the classical Dirac observables given above. Intro- 
duce operators Hq, f, h and to, which have, in particular, the commutators 



[h, Q = ih, [H ,h] = [f,h] = [Ho , t ] = [/, to] = 0. 
Then it is straightforward to verify that the operator 



F(t) := cxp 



-H (to + T) 



/exp 



-jH {t +T) 



commutes with the constraint, [F(t),Ho + h] = 0, and obeys the Heisenberg equation of motion 

d T F(r) = ^[H ,F(r)]. 



(3.28) 



(3.29) 



(3.30) 



In line with previous notation, we write F :— F(0). 

There is no guarantee a priori that the solutions to (|3.30[) will be simply related to the classical expressions 
(|3.24p - (|3.26p . We will insist on the commutation relations 

[v,p v ] = [(t>,p < t > }=ih, (3.31) 

with all other commutators amongst v,p vi 4> and p^ vanishing, and choose for the Hamiltonian the operator ordering 

1 



H = -6irGp v vp v + x p\. 



(3.32) 
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Then define 



V = +Pvv), 

which also implies p v = i(v~ 1 y + yv^ 1 ) and leads to the closed set of commutation relations 



We then indeed find 



[v, y] = ihv, [H , v] = ihl2-jrGy, [H , y] = -iHH . 

v(t) = v - 1 2wGyt - 67rGH t 2 
y(t) = y + H a t, 

from which V and Y can be defined by replacing t with to, and 

V(t) = V - VI-kGYt - %-kGHot 2 
Y(t) = Y + Hot. 



(3.33) 



(3.34) 



(3.35) 
(3.36) 



(3.37) 
(3.38) 



Since p^ commutes with Ho, we can also identify P*(t) = P$ = p$. The Dirac observable associated with <f> is 
more difficult to write down, as might be expected from the nonlinearity of the classical expression (|3.f 71) . However, 
this expression shows that e°"^^v(t), with a = \JYIttG, is quadratic in t, and it can be verified, by taking repeated 
commutators with H Q , that 



where 



m 



v + ato + -Ml 



(3.39) 



a = ap^-a \ y + — 



b = (a + a 2 y)v 1 a — a 2 Ho- 
Sadly, we have not found any useful closed-form expression for $(t). 



C. Physical Hilbert space 

Take a kinematical vector space of wavefunctions ^(v,p^,, h), in which operators v, p^ and h act by multiplication, 
while 



a* 

OV Opcf, Oil 



The Hamiltonian Ho is 



Ho^ = 



^a 2 d v vd v + \ 1 p 2 4) 



where a 2 = a 2 h 2 = YI-kGY? . and the constraint reads 



a 



The change of variables 



—vd v vd v + -p^ + vh 



2ip0 

5 

a 



%(v, P<t> ,h) = 0. 



2V2/i 



-v 



1/2 



a 



(3.40) 



(3.41) 



(3.42) 



(3.43) 
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converts this into a standard form of Bessel's equation 

(zd z zd z - v 2 + z 2 ) * = 0, (3.44) 

and we find 

2 

m = ^rz 2 ^ (3.45) 

Oil 

?)* = -j{zd z + l)ty (3.46) 



in* = -ihh,- 1 \^-zd z + hd h j * (3.47) 
ff * - A (^0, _ ^) (348) 

The constraint equation (|3.44[) is solved by 

*{v,P4,h)=il>fo,h)C v {z), (3.49) 

where C v {z) is some Bessel function, and ip(p<pi h) is arbitrary. (More generally, the solution is a linear combination 
of two independent Bessel functions of order v, but these two sectors are not mixed by any operator of interest, and 
we assume it is sufficient to keep just one.) Clearly, the physical Hilbert space is the space of functions ip(p,p,h), 
equipped with a suitable inner product. The operators fl3.45p - p.47p corresponding to 'partial observables', do not, of 
course, have any well-defined action in this space, but the Dirac observables do, and we easily verify that 

V^p = -~(a 2 d h hd h + h~%)ip (3.50) 

Yip = ihh 1/2 d h h 1/2 ip (3.51) 
H ^p = -hip. (3.52) 

The operator <3? is again harder to deal with. Note in particular that, although its conjugate momentum P$ is just 
multiplication by p^,, we cannot identify $ = ikd/dp^, since it does not commute with Hq. However, a well-defined 
action of the Dirac observable fl3.39[) can be found as follows. Acting on the Bessel function, we have 

e & k v (z) = e- 29 / d "C u (z) = C u - 2 (z). (3.53) 

Then, repeated use of the standard recurrence relations that connect C u ± n (z) and their derivatives suffices to verify 
that 

2 

e^V^v, h ) = ^ h - H v - 2 ) - - + Mdhhd h ] ^(v - 2, h), (3.54) 

where, as a shorthand, we use ip(v, h) to mean the same as ip(p<t>-i h) via (|3.43|) . 
Finally, make the change of variable 

h = ee Q " (3.55) 
where e is a positive constant that we will later choose to be close to (h). We find 

V = -le-^^+pJ) (3.56) 



a 

Ho = -ee Q " =: -h. (3.58) 

Choosing the inner product 

poo POO 

(^1,^2)=/ dp / di]e mi ^(p^,, rj)ip2 (PfaV), (3-59) 
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so 'Hphys = L 2 (M. 2 , dp^e^drj), these three operators and a fourth one derived from (|3 . 54[) are symmetric 6 . It is 
worth noting that this quantization scheme differs from the Wheeler-de-Witt scheme described in [13| not only by the 
addition of the clock variable, but also because the inner product is chosen to make the generator of displacements in 
r self-adjoint, rather than the generator of displacements in <j>. 

D. Viability of the test clock 

A wavefunction ip(p<f>,T}) in the physical Hilbert space specifies the state of a system consisting of two physical 
objects, the volume and the scalar field, represented by the pairs of Dirac observables (V, Y) and ($, P$) respectively, 
but the clock does not feature as an independent object; it has been eliminated by solving the constraint. As discussed 
in the previous section, we interpret this state of affairs by supposing that the clock is internal to an observer (or, less 
anthropomorphically, an observing apparatus) who is thereby debarred in principle from observing it. The remaining 
Dirac observables describe the volume and scalar field from the point of view of this observer, and evolve precisely as 
in standard quantum mechanics with respect to the 'heraclitian' time r. 

If this is to make sense, it must be possible to find states in which the clock serves to reveal the time evolution 
of the universe, while perturbing it to a negligible extent. One way of ensuring this would be to model the clock 
as an system whose energy is bounded. Thus, we might attempt to realize the Dirac observables (|3.50|) - (I3.52[) as 
self-adjoint operators on, say, W p hy S = £ 2 (R x [0, e], dp,/, dh). It is technically easier to retain our earlier choice 
'Hphys = L 2 (R 2 , dp^e^drj), on which the clock's energy h is positive (see (|3.55p ). and restrict attention to states in 
which this energy is small. It is not obvious that one of these choices is physically less reasonable than the other. 
Given the actual state of our world, the probability of finding the energy content of, say, a small alarm clock to be 
a significant fraction of the energy content of the visible universe is (presumably) extremely small. Depending on 
details of the clock's construction, a state in which that probability is appreciable might not be impossible in principle, 
though it would bear very little resemblance to the state that we actually have. 

In the model universe, we thus want to find states in which the probability of finding the clock's energy to be a 
significant fraction of the energy carried by the scalar field is very small; that is, we need both the expectation value 
e c := (h) and the variance (Sh) 2 :— ((h — e c ) 2 ) to be very small. Then the following issue arises. A simple-minded 
argument would suggest that, if 5h is vanishingly small, one has in effect an eigenstate of Hq, say Ho\e c ) = — £c\tc)- 
In that state, the expectation value of any commutator (e c | [F, Hq] |e c ) vanishes, and no time dependence should 
be apparent. This argument is in fact too simple-minded, because |e c ) is not a normalizable state in 'Hphys! what 
does happen is that the Dirac observables V(t) and Y(t) have very large variances when 5h is very small. Their 
probability distributions cannot be sharply peaked along some quasi-classical trajectory, and in that sense they are 
not well-defined time-dependent observable quantities. The important issue, then, is whether a useful compromise 
exists, such that the clock contributes, with high probability, only a negligible amount to the total energy of the 
universe, while the time-dependent observables remain reasonably well defined. 

In any state, the expectation value of the time-dependent volume operator p.37[) is 

V(t) := (V-(r)) = V + VHVt + ^e c r 2 , (3.60) 

where 

V:=(V) and H := -^(Y) (3.61) 

are the volume and Hubble parameter when the clock reads 0. (We emphasize that V(t) is the volume when an 
interval r of geometrical proper time has elapsed since that fiducial event: it cannnot be interpreted as the volume 
when the clock reads r.) The first two terms in V(t) follow the classical solution in the absence of the clock. When 
r is large enough, we have V(t) (a 2 e c j 'h 2 )T 2 , corresponding to a scale factor a(r) cx r 2 / 3 . In this situation, the 
energy density of the scalar field has become so dilute that the energy content of the fiducial cell is dominated by the 
clock, which in this model is equivalent to pressureless matter (see the remark following (|3.3jl ). At that point, the 
clock has ceased to qualify as a 'test' clock, but in a suitably chosen state, this will occur too far in the future to be 
a practical concern. 



We do not address the issue of self-adjointness; it is sufficient for our illustrative purposes that these operators have well denned matrix 
elements between Gaussian wavefunctions. 
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Consider a state of the form 



(3.62) 



where x(p<l>) 1S sufficiently sharply peaked at some value that p^, in (|3.56j) can be taken as just that value, which we 
call simply p$. The wavenumber q is real; an imaginary part corresponds to the constant e in (I3.55|) . The clock's 
energy is 



e c = (h) = ee 3 " 2 / 4A 



(3.63) 



and we define 



E A := -V- 



(3.64) 



which is the total scalar-field energy contained in the volume V. Then the expectation value of V turns out to be 



2 c 



1 



-A + q 



2 „2 



which can be rearranged to yield 



q 2 = 2VEs 



l + lL e -« 2 /2A 
Ed, 



-A, 



while the Hubble parameter defined by (|3.61[) is 



(3.65) 



(3.66) 



v = — iv\ = aq 

For the variances of the three operators (|3.56p - (|3.58[) , we find 



(3.67) 



t)~ = ( eQ2/2A - 1 



SY_ 
Y 

5h 



3a 2 /2A 



+C 



A 

p a 2 /2A 



A E,j, 
e c V e c 



1 + ^ e -a 2 /2A 
Ea, 



1 a 2 E t 



4> 



1 + J^ e -« 2 /2A 
2Ea, 



1_ 

64 \7y 



2 i r_^_ s 2 

e r .V 



If the clock energy e c is to contribute a negligible amount to the constraint, we need 



(3.68) 
(3.69) 

(3.70) 



Ea 



(3.71) 



and if the first term in (SV/V) 2 is to be small, we also need 

,.2 



cr 
2A <<L 



With these approximations, we get 



(3.72) 



SV\ a 2 A E& 1 / A 



V 



2A e c V 



(3.73) 



and if this is to be small, we finally need 
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If these conditions are met, then q 2 w 2VE^ and (5Y/Y) 2 w (A/4e c V)(e c /E^) is small. Also, the Hubble parameter 
given by p.67[) becomes 



f— ) 



2 2^=^ 



8ttG 



3 V 



(3.75) 



which is the classical Friedmann equation. 

With the approximations (|3.71|) . (|3.72p and ()3.74|) . the variance of the time-dependent volume operator gives 



f6V(n 



V V 



a 

2A 



1 - it (3Mr)! 



A 

e c V e c 



1 



2E, 



■(3-Hr) 



(3.76) 



where, on the LHS, V in the denominator is ( V), not V(r). As expected, this leads to the same physics as the version 
of the Wheeler-de-Witt equation studied in 13]. Namely, the singularity occurs when 3"Hr « — 1, at which point 
SV(t) is substantially unchanged from its value at r = 0. 

The compromise needed for the 'test clock' to make sense consists in satisfying the two inequalities 



« 2 * e 2 c V 

T« A «tr 



(3.77) 



simultaneously. For an essentially classical state at r = 0, this is indeed possible — and by a huge margin, as might be 
expected. In SI units, we have, first of all 



12ttG n /i 2 



3.1 x l(T 94 Jm 3 . 



(3.78) 



Although our universe is not dominated by the energy of a massless scalar field, for illustrative purposes we take 
— E^/V ~ 10~ 9 J m -3 (roughly the observed current energy density) and a clock with energy e c = 1 kg.c 2 ~ 10 17 J. 
Then 



= fc. „ 10«Jm 3 . 

E^ 



(3.79) 



Evidently, there should be no problem finding a A in the range (|3.77p . 

The time r c at which the scalar field energy is so diluted that the clock begins to perturb the evolution significantly 
is given by p.60p as 



-IT—SUV = — T , 



(3.80) 



where tq = (3H)' 1 is the age of the universe at the fiducial time r = 0. If V is the size of the observable universe, 
say V — 10 80 m 3 , we get t c ~ 10 54 to, so the last term in p.60p can safely be ignored for practical purposes, but that 
would clearly be true for much smaller volumes too. 

If A/e c V <C e c / E^, then the last term in (|3.73[) is negligible, and the sum of the remaining terms is minimized by 



A = 



la 2 e 2 V 



2E dj 



(3.81) 



namely the geometric mean of the two extreme values in (|3.77p . It is possibly interesting to consider this compromise 
in the context of the Robertson-Schrodinger uncertainty relation [36| 



(SV) 2 (SY) 2 >- ([V,Y]) 



({V-V,Y-Y}) 



(3.82) 



With the three approximations (|3.7ip . (|3 . 72[) and (|3.74l) . this inequality is saturated, both sides being given approxi- 
mately by 



-V 



1 



2A A E^ 

a 2 e r V e c 



(3.83) 



and the value (|3. 811) of A makes the commutator and anticommutator terms equal. 
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IV. LOOP QUANTUM COSMOLOGY WITH A TEST CLOCK 
A. Loop classical cosmology 

The cosmological model treated in arises, in effect, from quantization of the classical Hamiltonian 

H = Hgr av (v,a) + H ma ,ttei(v,P4,) (4.1) 
a 2 B 

-ffgrav = — ^— ^sm 2 (2er) [sincrjcoscr, — coscr{sincr, \v\}] (4.2) 

^matter = B^ft. (4.3) 

In terms of the volume v and its conjugate momentum p v dehned above, the new variables are given by 

^2BsgnO>, a = -^^p v (4.4) 

where B = /3h 3 / 2 a~ 3 and the numerical constant (3 = 2 1 / 2 3 5 / 4 7~ 3 / 2 is approximately equal to 48.24, when the 

Barbero-Immirzi parameter 7 is given the value 0.2375 . . . obtained from calculations of black-hole entropy in loop 
quantum gravity [33, [38j . Their Poisson bracket is 

{a,v} = tT 1 . (4.5) 



The motivation for this Hamiltonian is explained in detail in |12Ml5j (see also [39j]); here, we note only that the sign 
of the volume variable v corresponds to the orientation of a physical co-triad relative to a fiducial one, and that Ho 
reduces to the original form (|2.5|) in the limit h —> 0. 

By using the basic Poisson bracket (|4.5I) to evaluate those in (|4.2j) . we arrive at the Hamiltonian 

#lcc = -^H sin 2 (2a) + B\v\~ l p% (4.6) 

which defines the classical theory that we will refer to as 'loop classical cosmology'. It is straightforward to obtain 
the classical equations of motion and the pair of solutions 



v{t) = ± 

a 



1 a 4 B 2 2 
1 + 



1/2 



(4.7) 



valid on the constraint surface -ffLCC — 0j with the boundary condition that the minimum volume occurs at t = 0. 
The maximum density 

_ 2i? 2 P 2 _ q?B 2 V3 aR . 

^ ~ ^min " 2 ~ 167T 2 7 3 G 2 ^ 1 j 

is the same as that found in [l4l | for the quantum theory, while the characteristic time H/a 2 B is of the order of the 
Planck time. 

B. Dirac observables 

'Polymer' quantization schemes of the kind used in [l3l . fl4| promote the Poisson bracket (|4.5[) to the commutation 
relation 

[M = -6, (4-9) 

where t) is the holonomy \) = c 1CT , and take the operators f) and v to act in a Hilbert space of functions ^{v) : which 
have support on a countable subset of the real line, 

v^{v) = v^!{v), \)^>{v) = m{v + 1). (4.10) 

The gravity + matter system has the Hamiltonian 

Ho = ^C + B\v\-H% (4.11) 
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where, with a suitable ordering of non-commuting operators in the quantized version of (14.21) . the action of C is 

6^{v) = \v + 2|#(i/ + 4) + \v - 2|*(z/ - 4) - (\v + 2\ + \v- 2|)*(i/). (4.12) 

As before, we want to add to this system a clock, with Hamiltonian h which, along with its conjugate variable to, 
commutes with v and f). Again, the total constraint is Hq + h = 0, and we would like to introduce time-dependent 
families of Dirac observables (|3.29[) that commute with this constraint. However, a direct construction of these 
observables similar to that described in section IIIIBI is likely to be feasible only when their time dependence is 
polynomial (or perhaps has some other simple form). Here we adopt the following, somewhat less direct strategy. 

We consider an auxiliary Hilbert space H aux = Hq ® "Hciock where the Hilbert space associated with the clock is, 
say, % c iock = L 2 (R, dh). The Hilbert space W is a space of functions ^(v.p^), on which we again take p^ to act by 
multiplication, and on which we need the Hamiltonian Hq to be self-adjoint. This space can itself be decomposed as 
Wo = "Hgrav <8> Wrj), with = L 2 (R, dptf,). In order for H to be self-adjoint on Ho, we need both \u\ and C to be 
self-adjoint on Hgra-v- They must, in particular, be symmetric, and this requires the inner product 

(*i,*aW =X]*i( z/ )*2(^), 

V 

but again we will not attempt to establish their self-adjointness. This inner product is different from the one used in 
[mEl; namely (^i, ^2) = J2i> \ V Y& i{ v )^ i{v) , which is needed for self-adjointness of the operator = — C\P\ whose 
positive square root is the generator of displacements in the scalar field <f>, regarded as an internal 'time'. 

For any given value of p^, the Hamiltonian H Q (p^), obtained by replacing p^, in (|4. 1 1[) with its eigenvalue p$ is a 
symmetric operator on H STav . We assume (without attempting a rigorous proof) that it has a set of 8- normalized 
eigenfunctions ^ E (v, p^) 

H (pt)$ E (v, P $) = E<S> E {v,p$), ($ B , SuOgr-v = S(E - E'), (4.13) 
and that any function in H anyi can be expressed as 

^(u,p^h) = J dE$ E (v,p ! p)?p E (h,p l p), (4.14) 

with 



Mh,P*) = (*b(p«), *(P*, fe)) grav ■ (4-15) 

Suppose that the action of an operator / on Hq (equivalently, an operator on H aU x that acts as the identity on 'Hciock) 
can be specified by a kernel f(E, P( j,; E' \p'^). That is, 

fat{vM,h) = J dE^ E ( V , P4> )^\h, P4> ), (4.16) 

with 

i> [ i\h, V<t> ) = JdE'J dp^f(E, P4> ;E',p'^ E ,(h,p^. (4.17) 
The action of the corresponding Dirac observable F(t) defined in (|3.29[) is easily found to be 

^ {T)) {K P4> ) = !dE> ( dp'^ E - E ^/ h f(E,p^E\p'M E ,(h + E-E',p'). (4.18) 



If ty{v,p$, h) is a solution of the constraint equation, then the expansion coefficient in (I4.14p has the form 

iP E (h,P4)=tP{E, P4 ,)5(E + h), (4.19) 

the function ip(E,p^) specifying a particular solution. Because F(t) is a Dirac observable, its action (|4.18[) on one 
solution of the constraint equation yields another solution, specified by the function 

^ F( ^(E,p^) = IdE' ( dp^fiE^E'^e^-^'^E 1 ^). (4.20) 
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Evidently, the physical Hilbert space "H p hys is a space of functions ip(E,p^,), equipped with an inner product chosen so 
as to confer self-adjointness on some class of operators, and the action of the Dirac observables on "H p hy S is specified 
by (|4.20|) . If the kernel can be expressed in terms of a differential operator T as 

f{E,p^E\p^)=T{E'^,d E ,,d^)8{E-E')8{p 4> -p' <j> ), (4.21) 

then V> (F(r)) (£,£v) = e^/^Fe-^/V^,^), and we recover an algebra of differential operators analogous to (|3.50[1 - 
(I332)) . with h = -E. 

If / is constructed from i>, f) and p^, but does not contain cj), it can be construed as an operator on Hgrav that 
depends parametrically onp^. Its kernel has the form f(E^p^;E',p'^) = f{E 1 E' 1 p t j,)5{p l f ) — p'^), with 

f(E,E\ P<t> ) = ($E&),f*E>to))w ( 4 - 22 ) 
and we will use this expression to estimate the volume operator V(t). 



C. Eigenfunctions of H' {p^,) 

As described in [3, HH, the constraint equation Hq^ = that applies in the absence of the clock can be solved 
exactly if we restrict H gra v to functions ^{v) that have support only at v — An, n <E Z; the operator C in (|4,12|l 
clearly has a well-defined action on this restricted space. We will also restrict attention to this sector, but are able to 
obtain only approximate solutions to the eigenvalue equation (|4.13l) for nonzero eigenvalues E. 

With the definition 

$eO,2>4>) = -^-XE(n,Pt) = \n\xE{n,P4,) (4.23) 
the eigenvalue equation we wish to solve is 

^ [(n + l)(2n + l)x E (n + 1) + (n - l)(2n - l)x E (n - 1) - <in 2 x E (n)] + -^ XE {n) = E\n\x E (n), (4.24) 
where we suppress the argument p$ of \ E , and on taking the Fourier transform 

X E (n) = - / dke- i2nk XE (k), XE (k) = V e i2nk XE (n) (4.25) 
71 Jo 

J n=—oo 

this becomes 

B (a 2 smkd k smkd k +p 2 ) XE (k) = E \-2\d h \ XE (k). (4.26) 
Finally, the change of variable sin A: = 1/ cosh a; yields (with an obvious economy of notation) 

(d 2 x + \- 2 ) X e(x) =E\-2icoshxd x \xE(x), (4.27) 
where we have defined A = ajp^ and E = E / a 2 B. With this parametrization, the inner product in (|4.13l) is 

1 f°° 

($E,$E')grav = IT I dx COsh xd xXE\x)dxXE> (x). (4.28) 
47T J_ oc 

In the form (|4.27p . the eigenvalue equation is trivial when E = 0, but hard to solve when E is nonzero. We will 

obtain an approximate solution by means of an expansion in powers of A, but our analysis will be largely devoid of 
rigour. In practical terms, the approximation is likely to be quite good; for example, we have A ~ 10~ 123 for the state 
considered in section UlI Dl Consider the ansatz 

Xe(x) = -iAWexp [iA _1 a; + if E (x) + Xg E (x)] , (4.29) 

where f E {x) and g E (x) are real functions expressible as power series in A, 

f E (x) = f (x) + \f 1 (x) + ..., g E (x) = go(x) + X 9l (x) + . . . . (4.30) 
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The prefactor — iA is inserted for later convenience and Af is a normalization constant. We wish, of course, to determine 
the functions f E (x) and g E (x) by substituting this ansatz into (I4.27[) , but the absolute value of the volume operator 
v = — 2icosha;9 x presents a difficulty. It is not hard to show that this operator is positive on the space of functions 
whose Fourier transforms J dxe~ lu)x x{x) have support only for uj > and vice versa, and that these two spaces are 

orthogonal. Clearly, the function e lA x lies in the positive-volume space, and one might hope that the same is true of 
the trial function (|4.29[) . at least within the expansion in powers of A, but we are not able to prove this. We proceed 
by removing the absolute value symbol in (|4.27[) . and verifying a posteriori the positivity of v on the space spanned 
by the resulting set of approximate eigenfunctions. 
With this simplification, we find 

f E (x) = -E$mhx + ^E 2 (cosYvxsmhx + x) + 0(A 2 ) (4.31) 
g E (x) = i^coshz- ^£ 2 cosh 2 x + 0(A 2 ). (4.32) 

To make systematic use of this expansion, we construct the polar representation of the function 

ts(x) := d xXE {x) = p E (z)e iMx) , (4.33) 

with the result 

A ~ 



p E (x) = Af 



1 - ^£:cosha; + 0(A 2 ) 



(4.34) 



8 E (x) = A 1 x - Esinhx + ^^(coshxsinhx + x) + 0(A 2 ). (4.35) 

An important check on the consistency of our procedure is now to obtain the inner product 

1 f°° 1 f°° 

($b, $ E 0grav = -r / dxcoBhxZ E (x)ZE'(x) = — dx cosh x p E (x)p E > (x)e^ E <»>I . (4.36) 

4?r J_ oc 47T J_ oa 

Self-adjointness of Hq requires that this be proportional to 8{E — E'), which in turn requires 

(E - E')- 1 d x [e E (x) - e E ,(x)\ occoshxp E (x)p E ,(x), (4.37) 

in order that the final integral in (14.36[) reduce to J dse l( - E ~ E ' s . This can be checked order by order in A, and is readily 

seen to hold at the order of our explicit calculations. We find that $s')grav = S(E — E'), if the normalization 
constant in (|4~29]) is Af = {2/a 2 B) 1 / 2 . 

Matrix elements of the volume operator are given by 

1 f°° 

($ B , i>$e' ) gra v = ~r / dx cosh x £ E (x)^) (x) (4.38) 

47T 



with 



To lowest order, we find 



and thus 



4*°(a0 = dx [tX E (x)} = ~2id x [coBhxZ E (x)] . (4.39) 



^ { E \x) = 2X- 1 coakxt B (x) = 2X- 1 (l - a i B 2 d 2 E ) 1 ' 2 $ E (x) (4.40) 



($ B , z>$ £ 0grav = 2A- 1 (1 - a A B 2 d 2 E ,) 1/2 S(E - E'). (4.41) 



It follows, at this order of approximation, that v is positive, (^f, ^^ , ) gra v > on the space spanned by these eigen- 
functions; obviously, by starting with x E ( x ) = exp(— iA _1 .T + . . .) we would find a complementary space on which 
v is negative. The form of the series (|4.34l) and (|4.35l) suggests that a similar approximation might be obtained by 
expanding in the eigenvalue E, which we also want to be small, since it is minus the energy of a test clock. The 
leading terms of such an expansion do agree with those obtained above, but we have not been able to develop it in a 
wholly consistent manner. 
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D. Time evolution of the volume 



The matrix element (|4.40l) together with (|4.20[) - (|4.22[) yields the expectation value of the Dirac observable associated 
with v as 



<%)) 



1/2 



(4.42) 



where we observe an obvious correspondence with the classical solution (|4.7[) . 

We would again like to assess the viability of a test clock, by restricting the clock's energy — E to a small range of 
values around a mean value e c , which is itself small. As in section Till Dl we take the dispersion in p^ to be negligible 
and use the Gaussian wavefunction 



7T 



1/4 



-A{E+e c ) 2 /2 



(4.43) 



This state differs in detail from (|3.62[) ; roughly speaking, the two variances are related by A = A/a 2 e 2 . Taking into 
account the rescaling v = (2B)~ 1 \v\, we find for the Dirac observable V(t) associated with the volume v ~ a 3 

1/2 



V(r) := (V(r)) = v 



1 



dse 



1 



Ahs 



(4.44) 



where vo = p^/aB is the minimum ('bounce') volume attained in the classical solution, and rp = h/a 2 B w 6.8 x 10 45 s 
is of the order of the Planck time. In this state, the bounce occurs at r = 0, but including a factor Q- 1T o E / h m the 
wavefunction (|4.43j) would displace it to r = —tq. We also need 



T > = Vn 



Ahs 



Ah 2 



(4.45) 



If the clock's energy is to be restricted to a small range near e c , we need A ^ e~ 2 , or A > a 2 , which reproduces 
the first inequality in (I3.77[) . 

Now consider a late time, r ^ rp, at which it ought to be possible to set up a quasi-classical state. If Ah 2 /Vp is 
not too large (and we shall soon see that this quantity should be small), the volume can be approximated as 



V(r) = 251 

TP 



1 + ^ + 01- 



(4.46) 



This reproduces the second term of (|3.60[) , linear in r, but misses the last term proportional to r 2 which, as discussed 
earlier, is important at very late times when the clock energy dominates that of the scalar field. One would expect 
the large- volume evolution to be insensitive to the quantization scheme, and it is likely that the discrepancy indicates 
a failure of the approximations (I4.34[) and (|4.35[) for large values of E. With the approximation (|4.46l) . we estimate 
the dispersion in the volume at late times as 
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v 2 Ah 2 AE^ 



(4.47) 



and requiring this to be small reproduces the second inequality in (|3.77|) . At late times, therefore, the criterion for the 
notion of a 'test clock' to be viable is the same in the LQC and Wheeler-de-Witt quantization schemes, as it should 
be. 

At the 'bounce' volume, which occurs at r = 0, we might expect this criterion to be modified. Assuming that 
Ah 2 /rp is small, we obtain the approximation 



V(0) = v 



1 Ah 2 



_3_ Ah 2 
32 



(4.48) 



and 



18 



In order to have this dispersion small, and to have the clock's energy close to e c , we need 

e~ 2 « A « e p 2 , (4.50) 

where ep = h/rp ~ 1.6 x 10 10 J is about 10 times the Planck energy. Evidently, independently of the quantum 
state, the clock's energy e c must be much larger than ep. This condition might be regarded as satisfied by the 1kg 
clock considered in section IIIID| but not by a huge margin. Whether e c is much smaller than the scalar- field energy 
depends, of course on the particular state considered. For the quasi-classical state of a region of size 10 80 m 3 considered 
in section IPlI D[ the constant = 2V 2 p^ is about 10 151 Jm 3 . At its 'bounce' volume, vq, this region contains a scalar 
field energy = p 2 /2u = e P p t j ) /2a « 10 122 ep, so the condition e P ^ e c E^ is not hard to satisfy. 

V. DISCUSSION 

In section UH we identified two features of a widely adopted relational approach to time evolution in generally- 
covariant systems, which we believe to be shortcomings of this approach. In the context of simplified time- 
reparametrization-invariant models such as the cosmological model studied in this paper, evolution with respect 
to an 'internal time' is described by operators of the form (|2.8[) . These operators depend parametrically on some 
variable (in this case 4>o) associated with a dynamical variable (in this case a scalar field) which is regarded as a 
physical clock. The shortcomings (in our view) of this account are two-fold: (i) it provides no 'heraclitian' time that 
would serve to make sense of, for example, the everyday experience that a well-constructed clock reads '10s' about 
ten seconds after it read '0'; (ii) a parameter such as (f>o cannot be interpreted as an observed value of the chosen 
physical clock variable, because no dynamical quantity exists that might be observed to have this value, either as a 
function on the reduced classical phase space, or as an operator on the physical Hilbert space of the quantized theory. 

In sections IIIII and IIV1 we described, within two standard quantization schemes, a possible means of alleviating 
these difficulties by introducing a 'test clock' that is to be regarded as internal to some specific observer (or observing 
apparatus). In this way, we could construct classical Dirac observables, or 'evolving constants of the motion' (|3.9p and 
their quantum-mechanical counterparts (|3.29[) . These gauge- invariant quantities evolve with respect to a heraclitian 
time r, their evolution being governed exactly by the usual Hamilton or Heisenberg equations of motion. Like the 
physical clocks in other relational schemes, the test clock is not represented by any independent Dirac observable. 
Intuitively, this makes sense, insofar as the clock is internal to the observer from whose point of view the time evolution 
is described, and thus inaccessible to that observer. Problem (ii) above does not arise, because r has the textbook 
interpretation of the geometrical proper time that elapses along this observer's worldline, and does not correspond to 
any reading of the clock. 

It is worth emphasizing that more is at stake than interpretation. In constructing the physical Hilbert space, one 
naturally wants time evolution to be generated by a self-adjoint operator. As noted earlier, different inner products, 
and to that extent different quantum theories, are needed to confer self-adjointness on the generator of evolution in the 
scalar field <f> or on the original Hamiltonian Hq, which generates evolution in r. It seems to us that self-adjointness of 
Hq is a more natural requirement. Any specific choice of an internal time variable is essentially an arbitrary matter, 
and it seems unnatural that the resulting quantum theory should depend on that arbitrary choice. 

If the notion of a 'test clock' is to make good sense, it should be possible to restrict the clock's energy to be a 
negligible fraction of the total. This places limits on the resolution with which the values of time-dependent observables 
might be determined, but we found in sections 1111 Dl and IIV Dl that reasonable compromises are possible. That is, one 
can take the clock's energy to be reasonably small, while leaving the time-dependent observables reasonably sharply 
defined. 

As noted in section 1111 Al our construction of time-dependent Dirac observables is similar, but not identical, to 
relational constructions adopted by several previous authors, and it is worth comparing the two approaches in some 
detail. For classical systems, the original idea of Rovelli has been developed in considerable generality by Dittrich 
[IS EH (see (42| - |44l | for a somewhat different perspective), and the quantum theory obtained by reduced-phase-space 
quantization is discussed in [45I l46j. In models of the kind discussed here, these constructions amount to replacing 
the function to (a;) defined by (I3.5|) with a function ti(x, p) defined by 

r(x(x,h)) = p, (5.1) 

and the Dirac observable F(t) defined in (|3.8|) with 

F 1 (p) = f(wMx, P ))- (5.2) 

It is again straightforward to show that {ti,h} = —1 and that, for each value of the parameter p, Fi(p) Poisson- 
commutes with the total constraint Hq + Ii. Classically, before the constraint is imposed, Fi(p) seems to be interpretable 



19 



as the value of / when the clock reading r(x) has the value p. The operator (|2.8I) is a quantum implementation of 
this Dirac observable for the case that / is the volume and r{4>) = 4> the scalar field. Clearly, this is a different 
construction from the one we have advocated, and we have argued that it does not have a satisfactory interpretation: 
before imposing the constraint, p seems to be an observed value of the clock reading r, but no corresponding observable 
R exists in the final, constrained theory. 

There is, however, a special case in which the net results (|3.8|) and (|5.2p of these two different constructions are 
algebraically indistinguishable; namely, when the solution of the equation of motion for r is linear in t, so that 
t\(x,p) = —r(x) + p. This is true for a special kind of clock, whose Hamiltonian h = p r is just the momentum 
conjugate to r. In that case, the total constraint has the 'deparametrized' form H = H a + p r , and the replacement 
p r = —ihdp yields a Schrodinger-like equation. It can be argued (see, for example, (U [23| ) that the Hamiltonian 
h = p r , which is unbounded below, cannot describe any physical clock. 7 However, Brown and Kuchaf [29j have 
modeled a dust-filled universe by using a collection of scalar fields, of which one, say T, when the classical equations of 
motion are satisfied, is linear in the proper time along particle worldlines. This model has been studied more recently 
by several authors from the point of view of the relational formalism [H, HD, 0, E3] (see also [3(1 HH). 

In particular, Amemiya and Koike [H[ investigate the quantum dynamics of a homogeneous universe whose matter 
content comprises Brown-Kuchaf dust together with classical radiation and a cosmological constant. Their Hamilto- 
nian constraint has the form Ho(a,p a ) +pt, where pt, the momentum conjugate to T, is essentially the total energy 
of the dust content of a compact universe. Quantization amounts to replacing px with —iHOt and p a with —\hd a 
to obtain the Schrodinger-like equation ikdr^iT, a) — Hq(cl, —ihda)^ (T , a) . More formally, the relational formalism 
yields a family of Dirac observables A(p) through the prescription (|5.2I) . with ti(T,p) = p — T. Reduced-phase- 
space quantization then promotes these to 'Heisenberg-picture' operators A(p) on the physical Hilbert space and the 
corresponding Schrodinger-picture wavefunction ^(p, A) is governed by 

ihd p ^(p, A) = H (A, -ihd A )^{p, A). (5.3) 

We retain the symbol p here to emphasize that this variable is a value assigned to the scalar field T, whose classical 
equation of motion happens to have the solution T(t) — T(0) + 1, whereas our variable r is, up to a choice of origin, 
the arc length t. Again, there is no operator on the physical Hilbert space to represent the dust, so although p is a 
value assigned to T, it cannot be regarded as the result of a measurement of T, though in this special case it can be 
loosely associated with proper time by appealing to the classical equation of motion. Of course, the constraint implies 
that the energy of the dust is given by — Hq, but this ought not to be the whole story: cosmologists expect to verify 
the Friedmann equation through independent measurements of the Hubble parameter and the distribution of matter. 
This is a matter of some practical importance if the Brown-Kuchaf fields are used to model the actual non-relativistic 
matter in our universe. In [34j], the 'Schrodinger' equation (|5.3[) is solved for conditions corresponding to an energy 
density of radiation much larger than is observed, and yet their universe remains substantially unaffected by this 
radiation at the earliest times, when classical cosmology would lead one to expect a radiation-dominated universe. 
We suspect that this reflects a rather large value of (—Hq) corresponding to a large density of dust, whose presence 
can be inferred only from the behaviour of the scale factor. Presumably, one could arrive at a closer approximation 
to our universe by considering a state in which (—Hq) is small enough, but that would not solve the interpretational 
problem that the dust is not directly observable. 

Alternatively, one can try to regard the Brown-Kuchaf 'dust' as being truly unobservable. This is the view taken 
in [33j], and in [4?} these authors speculate that this dust might be a candidate for dark matter, since it interacts only 
gravitationally. 8 On this view, one has to postulate that the universe contains, in addition to ordinary detectable 
matter, another species which is unobservable in principle, not merely because it does not interact, but because it 
does not feature in the physical phase space. This seems implausible to us, because the dust appears in the first 
instance on the same footing as any other matter, and it is only as a matter of technical convenience that one chooses 
to solve the constraints by eliminating the variables associated with this, rather than some other species. 



This objection is irrelevant to our construction, because our clock is not required to provide a linear measure of proper time. As 
implemented in section fill I the clock has a positive energy (see H3.55J and in section llVI its energy was allowed to be negative, with 
small probability, purely as a matter of technical convenience. The objection [2l| that measured readings of a quantum-mechanical 
clock with energy bounded below are not guaranteed to increase monotonically with time also does not apply, because our r does not 
represent any such measured values. 

In fact, it is claimed in these papers that the dust must have a negative energy density, in order to yield a physical Hamiltonian that 
is positive, but that this does not matter, because the dust is unobservable. Amemiya and Koike [34| accept the negative Hamiltonian 
implied by a positive energy density for their dust, on the grounds that the original gravitational Hamiltonian is in any case negative, 
and we think they are right to do so, at least in the context of their homogeneous cosmology. 
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By contrast, the view proposed in this paper is that a specific observer, say 0±, will account for time-dependent 
observations made along her own worldline by describing the physical phase space in terms of variables that do not 
include a clock internal to herself. It does seem plausible to us that she cannot, in principle, make measurements of the 
reading of this particular clock. If the model is sufficiently detailed, her description of the physical phase space may 
include the variables corresponding to a clock that is internal to a second observer, O2, whose readings she might, in 
principle, be able to measure (if a socially acceptable procedure could be agreed on) . To account for time-dependent 
observations made along his worldline, however, O2 will coordinatize the physical phase space using variables that do 
not include his own internal clock, though they would include the clock internal to O1. 9 

We have, of course, substantiated this view of time evolution only in the context of a greatly simplified cosmological 
model. Within this model, the time parameter t, which has an unambiguous classical interpretation as the proper time 
that elapses along an observer's worldline, survives quantization as a gauge-invariant c-number parameter, because the 
lapse function N(s) in (|2.4p is not a dynamical variable, and is not promoted to an operator in the quantum theory. 10 
However, by considering only a comoving observer, we have avoided dealing explicitly with the dynamical variables 
x A '(A) in (|3.1|) . which more generally are needed to specify precisely what is meant, in the final quantum theory, by 
'an observer's worldline'. To deal adequately with a genuinely localized observer in a more general spacetime is a 
much more challenging enterprise, which we plan to explore in future work. 
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